Multi-order interference is generally nonzero 
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It is demonstrated that the third-order interference, as obtained from explicit solutions of Maxwell's equations 
for realistic models of three-slit devices, including an idealized version of the three-slit device used in a recent 
three-slit experiment with light (U. Sinha et al.. Science 329, 418 (2010)), is generally nonzero. The hypothesis 
that the third-order interference should be zero is shown to be fatally flawed because it requires dropping the 
one-to-one correspondence between the symbols in the mathematical theory and the different experimental 
configurations. 
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I. INTRODUCTION 

According to the working hypothesis (WH) of Refs. lUl-Ilt], quantum interference between many different pathways is simply 
the sum of the effects from all pairs of pathways. In particular, application of the WH to a three-slit experiment yields ^ 

O- /(r,000) = |v/i(r) + v/2(r) + r3(r)|', (1) 

|. where Xj/j with j = 1,2,3 represents the amplitude of the wave emanating from the jth slit with the other two slits closed and r 

O denotes the position in space. Here and in the following, we denote the intensity of light recorded in a three-slit experiment by 

c/3 /(r, OOO), the triple O's indicating that all three slits are open. We will write /(r, COO) for the intensity of light recorded in the 

O experiment in which the first slit is closed, and so on. 

c/3 ' Assuming the WH to be correct, it follows that 

/(r,000) - |v/i(r) + v/2(r) + r3(r)P 

I r = |V/i(r) + v/2(r)P + |V^i(r) + VA3(r)|2 + |v/2(r) + r3(r)P 

(M ; -\m{r)\'-\V2{r)\'-\W3{r)\' 
>'. = /(r , OOC) + /(r, OCO) + /(r , COO) 

^! -/(r,OCC)-/(r,COC)-/(r,CCO). (2) 

I In other words, still assuming the WH to be correct, we must have 

A(r) = /(r,000)-/(r,OOC)-/(r,OCO)-/(r,COO) 
2 ! +/(r, OCC) + /(r, COC) + /(r, CCO) =0. (3) 



■ 
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According to Refs.lUS the identity Eq. Q follows from quantum theory and the assumption that the Born rule /(r) |1'(r) p 
holds. In a recent three-slit experiment with light [3], the seven contributions to A(r) were measured and taking into account the 
unceitainties intrinsic to these expeiiments, it was found that A(r) w 0. This finding was then taken as experimental evidence 
^ . that the Bom rule /(r) oc |*I'(r)p is not violated [3]. 

The purpose of the present paper is to draw attention to the fact that within Maxwell's theory or quantum theory, the premise 
that Eq. ([T]i (which implies Eq. (O) holds is generally false. By explicit solution of the Maxwell equations for several devices, 
including an idealized version of the three-slit device used in experiment [3J, we show that A(r) is generally nonzero. We also 
point out that the reasoning that leads to the WH and to the conclusion that A(r) = is fundamentally flawed because 



'Electronic address: 'h.a.de.raedt@rug.nl | 
^Electronic address: k.michielsen@fz-juel ich.de"| 
^Electronic address: k-hess@illinois.edu „ 



2 




FIG. 1: Amplitudes of the (left) and E. (right) components of the electric fields as obtained from a FDTD solution of Maxwell's equation 
for light incident on a metallic plate with three slits. The incident wave is monochromatic and has wavelength X = 500nm. The slits are X 
wide, their centres being separated by 3A. The index of refraction of the 4A-thick metallic plate (colored black) is 2.29 + 2.61/. In the FDTD 
simulations, the material (steel) is represented by a Drude model |4]. 
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FIG. 2: Left: Angular distribution of light 7(0; 000) transmitted by = 3 slits (see Fig.[T]( as obtained from the FDTD simulation (bullets) 
and Fraunhofer theory (solid line) 1(6, s = l,d = 2,N = 3), see Eq. l|6j, where s and d are the dimensionless slit width and slit separa- 
tion, respectively |5]. Right: Normalized difference £(6) = [/(0;OOO) -/(e;COO) -/(e;OCO) -/(0;OOC) +/(e;CCO) +/(0;COC) + 
/(0;OCC)]//(e = 0;OOO) as a function of 6. According to the WH of Refs. [IHl, this difference should be zero. 



it does not satisfy one of the basic criteria of a proper mathematical description of a collection of experiments, namely that 
there should be a one-to-one correspondence between the symbols in the mathematical theory and the different experimental 
configurations. 



II. SOLUTION OF MAXWELL'S EQUATION 

A conclusive test of the WH Eq. ^ is to simply solve the Maxwell equations for a three-slit device in which slits can be 
opened or closed (simulation results for the device employed in the experiment reported in RefjUare presented in Section lTVT) . 
For simplicity, we assume translational invaiiance in the direction along the long axis of the slits, effectively reducing the 
dimension of the computational problem by one. 

In Fig. [1] we show the stationary-state solution of the Maxwell equations, as obtained from a finite-difference time-domain 
(FDTD) simulation |4]. From the simulation data, we extract the angular distribution 7(0,000). Repeating these simulations 
with one and two of the slits closed, we obtain 7(0. COO) and so on. In all these simulations, the number of mesh points per 
wavelength A was taken to be 100 to ensure that the discretization errors of the electromagnetic (EM) fields and geometry are 
negligible. The simulation box is 75A x 40A large (coiTesponding to 30 01 1 501 giid points), terminated by UPML boundaries 
to suppress reflection from the boundaries [4]. The device is illuminated from the bottom (Fig.[T]i, using a cuiTent source that 
generates a monochromatic plane wave that propagates in the vertical direction. 

In Fig.|2tleft) we show a comparison between the angular distribution of the transmitted intensity 7(0;OOO) as obtained from 
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the FDTD simulation (bullets) and Fraunhofer theory (solid line). Plotting 

/(0;OOO) -/(0;COO) -/(0;OCO) -/(0;OOC) +/(0;CCO) 



^ ' /(0=O;OOO) 
/(0;COC)+/(0;OCC) 

(4) 



/(0 = O;OOO) 

as a function of (see Fig. |2lright)) shows, beyond any reasonable doubt, that the WH Eq. ([TJ of Refs. lU-dl] is in conflict with 
Maxwell's theory: E(0) takes values in the 0.5% range, much too large to be disposed of as numerical noise. Note that E(0) is 
obtained from data produced by seven different device configurations. 

Physically, the fact that E(0) ^ is related to the presence of wave amplitude in the vicinity of the surfaces of the scattering 
object (one, two, or three slit system), see for instance Fig. [itright). These amplitudes are very sensitive to changes in the 
geometry of the device, in particular to the presence or absence of a sharp edge. Although these amplitudes themselves do 
not significantly contribute to the transmitted light in the forward direction, it is well-known that their existence affects the 
transmission properties of the device as a whole (mItI- 



A. Illustrative example 

The essence of a wave theory is that the whole system is described by one, and only one, wave function. Decomposing this 
wave function in various parts that are solutions of other problems and/or to attach physical relevance to parts of the wave is 
a potential source for incorrect conclusions and paradoxes. Even for one-and-the-same problem, the idea to think in terms of 
waves made up of other waves can lead to nonsensical conclusions, such as that part of a light pulse can travel at a superluminal 
velocity. Of course, we may express the wave field as a superposition of a complete set of basis functions, e.g. by Fourier 
decomposition, and this may be very useful to actually solve the mathematical problem (to a good approximation). However 
such decompositions are primarily convenient mathematical tricks which, in view of the fact that in principle any complete set 
of basis functions could be used, should not be over-interpreted as being physically relevant 

The WH Eq. ([TJ takes these ideas substantially further by decomposing the wave amplitude in three parts, each part describing 
the same system (a single slit) located at a different position in space. It is then conjectured that the amplitude for the whole 
system (three slits) is just the sum of these three different amplitudes. 

Advocates of the "physical" motivation for this conjecture might appeal to Feynman's path integral formulation ifioll of wave 
mechanics to justify their picture but in fact, one can see immediately from Feynman's path integral formalism that the WH 
Eq. ([T]i is not valid. 

We use the expression for the propagator of the electron as given by Feynman and assume that the particle proceeds from a 
location a and time f„ on one side of the screen with slits labeled 1,2,3 to a location h where a measurement is taken at time t\, 
on the other side. We assume that there exists some time tc between ta and (as assumed by Feynman on p. 36, Ref.[T(|). The 
propagator for this process is denoted by Feynman as K{b.a) ifToll . If we include for clarity the times then we would have to 
write K{{b,tij), (fl,f„)). As pointed out by Feynman (Ref.|T^ p. 57) we have a connection of this propagator to the wave function 
i/A given by: 

\l/{b,tb)^K{{b,th),ia,ta)). (5) 

Feynman represented the propagator Khy a path integral that sums over all possible space-time paths to go from ato b with the 
end-point times as given above. If we have an infinitely extended screen in between a,b with only slit 1 open, then all paths can 
only proceed through this one slit. We denote the wave function that is calculated for a path leading through a particular point 
xi of the slit at time fvi by y/f. Similarly for slits 2 and 3 open only we have ™d 1/3 respectively and the corresponding K's 
are calculated with Feynman paths that only go through slits 2 or 3 respectively. 

Had we chosen all three slits open, then Feynman's formalism insists that pathways going through multiple slits matter in 
general. Therefore, we would have to include paths through multiple slits in the path integral representation of K and we would 
obtain a corresponding ^^[2^- Thus, Feynman's quantum mechanics with all three slits open does contain an infinity of paths that 
go through multiple slits resulting in ^'123- However, none of the wave functions \//[, or V3 ^nay contain any path through 
more than one slit because of the assumption that only one slit be open at a time. Therefore all the expressions involving these 
amplitudes do not contain multiple-slit path integrals and consequently do not contain all the paths that are required to compute 
'^123- illustrate the importance of "all" by solving the Maxwell equations for a minor variation of the three-slit experiment 
in which we block one slit. 

The geometry of the device that we consider is depicted in Fig. [3] together with the FDTD solution of the EM fields in the 
stationary state. We have taken the three-slit device and blocked the middle slit by filling half of this slit with material (the 
same as used for other parts of the three-slit device), once from the top. Fig. [3jtop), and once from the bottom. Fig. |3lmiddle). 
Comparing the FDTD solutions shown in Fig.|3ltop) and Fig.[3lmiddle), it is obvious to the eye that the wave amplitudes provide 



4 




FIG. 3: Top and middle: Amplitudes of the Ex (left) and £, (right) components of the electric fields as obtained from a FDTD solution of 
Maxwell's equation for light incident on a metallic plate with two slits and a hole between the two slits. The incident wave is monochromatic 
and has wavelength A = 500nm. The slits are X wide, their centres being separated by 3A. The index of refraction of the 4A -thick metallic 
plate (steel, colored black) is 2.29 + 2.61;. The holes are A wide and 2A deep. Bottom left: Angular distribution of light transmitted by the 
devices extracted from the FDTD simulation data. Solid line: Slit in the centre filled with material half-way from the top (see top panel); Red 
bullets: Slit in the centre filled with material half-way from the bottom (see center panel). On the scale used, the two angular distributions 
cannot be distinguished. Right: Normalized difference between the angular distributions of the device with the hole in the bottom (top panel) 
and the top (middle panel). Adopting the reasoning of Refs. Q,-i3J, this difference should be zero. 



no support for the idea that these systems can be described by a wave going through one slit and another wave going through the 
other slit. As in the case of the three slits, the angular distributions for the two cases look very similar (see Fig. |3jbottom,left)) 
but differ on the one-percent level. 



III. WHEN DOES THE WORKING HYPOTHESIS HOLD? 

It is of interest to scrutinize the situations for which the WH Eq. ([T]) is correct ifTTl [T2I1 . As pointed out earlier, in general, 
interference between many different pathways is not simply the sum of the effects from all pairs of pathways. To establish 
nontrivial conditions under which it truly is a pairwise sum, we discard experiments for which the WH trivially holds, that is 
we discard experiments that exactly probe the interference of three waves, such as the extended Mach-Zehnder interferometer 
experiment described in Ref. 2 and the class of statistical problems described by trichotomous variables considered in Ref. 13. 

Let us (1) neglect the vector character of EM waves and (2) assume that the diffraction of the three-slit system is described 
by Fraunhofer diffraction theory. Then, for normal incidence, the angular distribution of light intensity produced by diffraction 
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from slits is given by ^ 

( sinfA^TTt/sin 0) \ ^ / sinfTT^sin 0) \ ^ 

I{0,S,d,N = .\ . ' 1 ■ n ) ' (6) 

where s and d are the dimensionless slit width and slit separation expressed in units of the wavelength X, respectively. Therefore, 
we have 



A(0) = I{e,s,d,3)-2I{e,s,d,2)-I{e,s,2d,2) + 3I{e,s,d,l) 

= [(1 +2cos2fli/)2-8cos2flt/-4cos22fl(i + 3] ("^^^^ =0, (7) 



where a = tt sin 0. Thus, in the Fraunhofer regime the WH Eq. ([T} holds. 

It is not difficult to see that A(0) = is an accident rather than a general result by simply writing down the Maxwell curl 
equations Sill] 

e(r)^^ = VxH(r,0-/(r,0 

M(r)^?^ = VxE(r,f), (8) 

where the geometry of the device is accounted for by the permittivity e(r) and for simplicity, as is often done in optics Jsj], we 
may assume that the permeability (r) = 1 . 

Let us write e(r,000) for the permittivity of the three-slit geometry and E(r,f,000), H(r,f,000) for the corresponding 
solution of the Maxwell equations Eq. (O. The WH Eq. ([T]) asserts that there should be a relation between (E(r,f,000), 
H(r,f ,000)) and (E(r,f ,COO), H(r,f ,COO)), (E(r,f,OCO), H(r,f ,OCO)), (E(r,f ,OCC), H(r,f,OCC)) but this assertion 
is absurd: There is no theorem in Maxwell's theory that relates the solutions for the case e(r,000) to solutions for the cases 
e(r,COO)... e(r,OCC). The Maxwell equations are linear equations with respect to the EM fields but combining solutions for 
different e's is like adding apples and oranges. 

Of course, this general argument applies to the Schrodinger equation as well. For a particle moving in a potential, we have 



ifi^^ = (^-p" + V(r)j'P(r,0. (9) 

In essence, the WH Eq. ([T]) asserts that there is a relation between the solutions of four problems defined by the potential V (r) 
and three other potentials V; (r) for j — 1,2,3. More specifically, it asserts that 

^(r,f) = ^i(r,r)+^2(r,f) + ^3(r,f), (10) 

and 

y(r)1'(r,0 = y(r)^i(r,0+y(r)1'2(r,f)+V(r)^3(r,f) 

= yi(r)^i(r,f)+y2(r)^2(r,f)+y3(r)^3(r,0- (H) 

The authors could not think of a general physical situation that would result in Eq. (fTTT i. 



IV. DOES THE EXPERIMENT REPORTED IN REE J, HAVE MERIT? 

Although in general, the WH of Refs. fT-^ does not hold, the three-slit experiment reported in Ref. lHjlis not without merit, 
as we now show by discussing FDTD simulations of an idealization of the device used in the experiment |3|]. 

The geometry of this device is depicted in Fig. |4] (see also Ref. HI), together with the stationary state FDTD solution of the 
Maxwell equations. In the simulation, the device is illuminated from the bottom (Fig.[T]i, using a current source that generates a 
monochromatic plane wave that propagates in the vertical direction. The wavelength of the light, the dimension of the slits and 
their separation, blocking masks and material properties are taken from Ref. dH]. In view of the large (compared to wavelength) 
dimensions of the slits, to reduce the computational burden, we assume translational invariance in the direction along the long 
axis of the slits. This idealization of the real experiment does not affect the conclusions, on the contrary: It eliminates effects of 
the finite length of the slits. In all these simulations, the 81 mesh points per wavelength (A = 405 nm) were taken to ensure that 
the discretization errors of the EM fields and geometry are negligible. The simulation box of 820/xm x 120/xm (corresponding 



FIG. 4: Two-dimensional representation of the experiment reported in Ref. [3]. The three slits at the top are 30/im wide, their centres being 
separated by 100/xm. The blocking mask at the bottom can have one, two or three slits, each slit being bOflm wide with its centre aligned 
with one of the slits in the top plate |3]. In the example shown, the middle slit of the blocking mask is closed (corresponding to the case 
OCO). The separation between the top plate and blocking mask is 50flm. The index of refraction of the 25/xm-thick material (colored black) 
is 2.29 + 2.61i (index of refraction of iron at 405nm). The wavelength of the incident light is 405nm. Also shown are the amplitudes of the 
Ex (left) and (right) components of the electric fields as obtained from a FDTD solution of Maxwell's equation for a monochromatic light 
source (not shown) illuminating the blocking mask. Note that the Ex- and £;-components propagate in a very different manner. 




FIG. 5: Left: Angular distribution of light transmitted by the system shown in Fig. |4] for the cases in which all slits are open (000), 
one slit is closed (COO,OCO,OOC), and two slits are closed (CCO,COC,OCC), as obtained from FDTD simulations. Right: Normalized 
difference £(6) = [/(e;000) -/(e;COO) -/(e;OCO) -/(e;OOC) +/(0;CCO) +/(e;COC) +/(e;OCC)]//(e = 0;OOO) as a function 
of 6. According to the WH 1 1, 3], this difference should be zero. Note that the symmetric structure of Z(6) rules out that the signal is due to 
numerical noise. 



to 3 936 188 001 grid points) contains UPML layers to eliminate reflection from the boundaries [4]. Each calculation requires 
about 900GB of memory and took about 12 hours, using 8192 processors of the IBM BlueGene/P at the Jiilich Supercomputing 
Centre. 

Qualitatively, Fig.|4jleft) indicates that the x-component (Ex) of the EM-field propagates through the two layers of slits with 
very little diffraction from the top (=blocking) layer This is not the case for the z-component shown in Fig. Iljright). In this 
idealized simulation setup, the amplitude of the y-component of the EM-field is zero. 

In Fig.|5l;ieft) we present the results for the angular distribution of the seven cases (OOO, OOC, OCO, COO, OCC, COC, and 
CCO), extracted from seven FDTD simulations (curves with two slits closed overlap). From Fig. ISjright), it is clear that 
is not identically zero. However, is of the order of 10^^ and such small values of are unlikely to be observable in an 
experiment such as the one reported in Ref. 3^- Another way of phrasing this is to say that the experiment [3J is performed in a 
regime in which scalar Fraunhofer theory works well, as can be expected from the dimensions of the slits and slit separations 
of the device. As we have shown in Section [nil it is precisely in this limit that the WH holds. Therefore, assuming that the 
non-linearity of the detectors used in the experiment [3] is negligible, the experiment [3J suggests that the measured intensity is 
proportional to the square of the amplitudes of the wave field interacting with the detector. 



V. DISCUSSION 



A necessary condition for a mathematical model to give a logically consistent description of the experimental facts is that there 
is one-to-one correspondence between the symbols in the mathematical description and the actual experimental configurations. 
When applied to the three-slit experiment in which one slit or two slits may be closed, the argument that leads from Eq. ([T) to 
Eq. ^ is fundamentally flawed in that there is no such correspondence. 

If xj/j in Eq. ^ is to represent the amplitude of the wave emanating from the jth slit with all other slits closed, the WH should 
be written as 

/(r,000) = |v/i(r,OCC) + r2(r,COC) + v/3(r,CCO)p, (12) 
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that is, we should label the y/s such that there can be no doubt about the experiment that they describe. This notation establishes 
the necessary one-to-one correspondence between the mathematical description (the y/'s) of the particular experiment (labeled 
by OCC, etc.). Now, we have 

/(r,000) = |v/i(r,OCC) + v/2(r,COC)|2 + |v/i(r,OCC) + v/3(r,CCO)|2 
+ |VA2(r,COC) + VA3(r,CCO)|2 

-|V/i(r,OCC)|2-|v/2(r,CCO)|2-|v/3(r,CCO)p. (13) 

At this point, it is simply impossible to bring Eq. (fT3T l into the form Eq. (|2]i without making the assumption that 

V/(r,OOC) = i/i(r,OCC) + v/2(r,COC), 
VA(r,OCO) = v^i(r,OCC) + v/3(r,CCO), 

V/(r,COO) = i/2(r,COC) + v/3(r,CCO). (14) 

If we accept this assumption, we recover Eq. (|2]l. However, the assumption expressed by Eq. (fl4l i cannot be justified from 
general principles of quantum theory or Maxwell's theory: The only way to "justify" Eq. (fT4l i is to "forget" that the y/'s are 
labeled by th e ty pe of experiment (e.g. OCC) they describe. For a discussion of this point in the case of a two-slit experiment, 
see Refs. [Hill. 

In other words, accepting Eq. (fTTt destroys the one-to-one correspondence between the symbols in the mathematical theory 
and the different experimental configurations, opening the route to conclusions that cannot be derived from the theory proper 
Hence, if A(r) would indeed be different from zero for a three-slit experiment, one cannot conclude that Bom's rule does not 
strictly hold. However, if A(r) is found to be nonzero in the case of the Mach-Zehnder interferometer experiment described in 
Ref.^and non-linearity of the apparatus can be ruled out, the conclusion would be that Born's rule does not hold. 

VI. SUMMARY 

The results of this paper can be summarized as follows: 

1 . The third-order interference, as obtained from explicit solutions of Maxwell's equations for realistic models of three-slit 
devices, is nonzero. 

2. The hypothesis that the third-order interference Eq. (O should be zero is fatally flawed because it requires dropping 
the one-to-one correspondence between the symbols in the mathematical theory and the different experimental configura- 
tions. 

3. Having shown that the expression for the third-order interference Eq. (|3) cannot be derived from quantum theory or 
Maxwell theory of a three-slit experiment, it follows that any conclusion based on this expression has nothing to say about 
these theories (with the exceptions mentioned in Sectionllllli. 
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